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1.0  INTRODUCTION 

The  verification  of  the  gas  dynamic,  radiative,  and  contaminant  properties  of  high 
altitude  rocket  exhaust  plumes  requires  the  use  of  ground-based  simulation  facilities.  The 
accurate,  high  altitude  simulation  of  these  phenomena  is,  if  not  impossible,  often  quite 
difficult,  since  the  high  enthalpy  volume  flow  rate  for  a  multisecond,  steady-state  firing  can 
exceed  greatly  the  effective  pumping  speed  of  the  typical  cryopumped  test  chamber.  An 
obvious  solution  to  this  problem  is  the  use  of  short-duration,  pulsed  engine  firings.  For  this 
mode  of  operation  the  combined  advantages  of  the  nonsteady-state  behavior  of  the 
interaction  of  the  exhaust  plume  and  the  proportionately  decreased  time-integrated  volume 
flow  rate  into  the  chamber  serve  to  maintain  low  chamber  pressures  during  the  engine  firing 
and  to  enable  rapid  chamber  pumpdown  between  firings.  Consequently,  simulated  high 
altitude  testing  capabilities  are  achieved  without  the  necessity  of  using  a  significantly  larger 
cryopumped  chamber. 

Clearly,  to  achieve  this  capability  of  pulsed  engine  testing  the  required  diagnostics  must 
possess  adequate  speed  and  sensitivity  to  yield  the  necessary  spatial  and  temporal 
resolutions,  and,  further,  the  nonsteady  behavior  of  the  chamber  performance  must  be 
understood  to  define  adequately  the  simulated  altitude-versus-time  characteristics.  Both 
areas  are  under  active  study  at  the  Arnold  Engineering  Development  Center  (AEDC). 
However,  an  additional  obstacle  for  attaining  this  capability  is  ensuring  that  the  pulsed 
engine  exhaust  plume  of  a  liquid  fuel  engine  is  a  faithful  representation  of  the  multisecond 
firing  engine  exhaust;  i.e.,  steady-state  characteristics  must  be  achieved.  Order-of- 
magnitude  calculations  and  experimental  data  of  other  work  can  be  provided  for  typical 
engines  to  show  that  steady  state  is  definitely  achieved  for  development  of  the  nozzle  flow, 
gaseous  species  combustion,  and  radiative  emission  properties.  However,  central  to  the 
feasibility  of  pulsed  engine  testing  for  contamination  studies  is  the  question  of  whether 
steady-state  formation  of  particulates  within  the  exhaust  plume  is  achieved  in  short  firing 
times.  Specifically,  for  particulates  formed  by  flow-field  condensation  processes  this 
question  can  be  addressed  analytically.  However,  the  formation  of  particulates  by  the 
incomplete  combustion  of  fuel  droplets  is  much  more  difficult  to  predict  theoretically,  and 
recourse  must  be  made  to  experimental  studies.  The  goal  of  this  work  was  to  determine 
whether  steady-state  particulate  formation  existed  for  short-duration  (20-  to  100-msec) 
firings  of  a  5-lbf-thrust  liquid  bipropellant  engine,  the  exhaust  plume  properties  of  which 
have  been  recently  studied  (Refs.  1  and  2)  in  the  VKF  Aerospace  Chamber  (10V). 
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2.0  EXPERIMENTAL  CONFIGURATION  AND  THEORY 

During  this  just-mentioned  study,  as  described  in  Ref.  2,  laser  scattering  was  used  as  a 
flow  visualization  technique.  The  Rayleigh/Mie  radiation  scattered  by  the  flow  field  from 
incident  pulsed  ruby  laser  beam  was  detected  by  an  image  intensifier-vidicon  system.  The 
laser  beam  was  injected  into  the  flow  field  perpendicular  to  the  flow  axis  at  a  distance  of  1 . 1 
in.  (2.8  cm)  from  the  nozzle  exit  and,  using  the  intensifier,  observations  in  the  radial 
direction  of  ±  0.735  in.  (±1.88  cm)  relative  to  the  centerline  were  performed  with  a  spatial 
resolution  of  0.1  in.  (0.25  cm).  The  image  intensifier  viewed  the  plume  at  90  deg  with  respect 
to  both  the  illuminating  laser  beam  and  the  plume  axis,  and  the  incident  laser  polarization 
axis  was  oriented  very  nearly  antiparallel  to  the  scattering  plane  (plane  of  the  incident  and 
scattered  light  propagation  vectors). 

Data  were  acquired  with  a  nominal  time  resolution  of  1  msec.  For  the  majority  of  the 
data  acquired  the  laser  scattering  signal  was  unpolarized.  However,  during  the  investigation 
of  an  engine  exhaust  combustion  pulse  of  a  nominal  20-msec  duration,  the  degree  of 
polarization  of  the  scattered  radiation  was  determined  at  the  times  t  —  17  msec  and  60  msec, 
the  latter  point  occurring  in  the  post-shutdown  phase  of  the  burn.  (The  times  quoted  in  this 
work,  as  in  Ref.  2,  are  relative  to  the  leading  edge  of  the  electrical  pulse  which  opens  the 
propellant  valves.) 

Using  the  experimental  apparatus  described  in  Ref.  2  and  a  calibration  procedure  which 
used  the  Rayleigh  scattering  from  N2  at  a  known  temperature  and  pressure,  the  unpolarized 
Rayleigh/Mie  scattered  radiation  was  used  to  provide  spatially  and  temporally  resolved 
values  of  the  product  of  number  density,  n,  and  scattering  cross  section,  a(x),  corresponding 
to  the  particle  size  distribution  function  (PSDF)  of  the  effluent.  The  average  size  parameter, 
(x),  is  defined  in  terms  of  the  mean  diameter,  D,  and  laser  wavelength,  X,  to  be 


X  !_  77  1}  A 

Since  the  purpose  of  the  laser  scattering  measurements  of  the  exhaust  plume  was  flow 
visualization  rather  than  particulate  characterization,  insufficient  data  exist  to  determine 
independent  values  of  n  and  a,  which  depend  on  the  particulate  mean  size  parameter  and 
index  of  refraction,  r\.  Consequently,  no  must  serve  as  a  measure  of  the  time  dependence  of 
the  particulate  properties;  i.e.,  if  no  is  found  to  be  constant  in  time  it  will  be  assumed  that 
both  n  and  a  (and,  therefore,  D)  are  all  constant.  It  is  recognized  that  it  is  possible,  in 
principle,  that  a  constant  na  product  can  be  obtained  for  time-varying  values  of  n  and  a,  but 
such  a  coincidence  is  assumed  not  to  be  the  case.  Moreover, for  particle  size  distribution 
functions  of  comparatively  narrow  width,  it  can  be  shown  that  for  x  values  on  the  order  of  2 
and  less,  with  approximately  a  factor  of  two  (2)  uncertainty, 
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<7  =  D3 

In  the  absence  of  information  on  the  shape  of  the  PSDF,  we  assume  that  D3  =  (D)3  and  that 
D  will  be  obtained  from  the  analysis  of  the  polarized  scattered  light.  Thus, 

no  -  pp 

where  qp  is  the  particulate  mass  density  within  the  flow  field.  Therefore,  one  can  regard  the 
measured  values  of  na  as  a  measure  of  the  mass  density  of  particulates.  It  is  seen  that  a  time- 
independent  behavior  of  n<r  implies  a  constant  value  of  gp,  from  which  time-independent 
values  of  n  and  D  are  assumed. 


The  polarization-specific  data  acquired  at  17  msec  in  one  set  of  measurements  were  used 
to  determine  (Refs.  2  and  3)  the  average  size  parameter,  x.  Assuming  x  to  be  constant 
throughout  the  engine  combustion  period,  local  values  of  n  were  determined,  which,  when 
combined  with  estimates  of  the  flow  speed,  v^,,  yielded  estimates  of  the  particulate  flux 
density,  Jp,  where 

<JP  =  nv~ 


3.0  RESULTS  AND  DISCUSSION 
3.1  LASER  SCATTERING  RESULTS 

Shown  in  Figs.  1  through  7  are  the  variations  with  engine  firing  time  of  the  nu  product 
for  the  radial,  or  vertical,  region  0.735  in.  (-1.88  cm)  sys  0  in.,  and  it  is  seen  that  for  all 
off-axis  locations  the  product  of  number  density  and  average  cross  section  reaches  a  steady- 
state  value  in  less  than  30  msec.  Further,  for  radial  positions  |y[  <  0.625  in.  (1.57  cm), 
steady  state  is  attained  in  less  than  20  msec.  Similar  results  were  found  for  ail  other  cases 
studied.  Insofar  as  this  project  is  concerned,  the  project  goal  has  been  satisfied  at  this  point 
in  that  the  particulate  properties  as  manifested  by  the  n a  product  have  been  shown  to  reach  a 
steady-state  condition  in  times  on  the  order  of  20  to  30  msec.  However,  it  is  of  interest  with 
regard  to  the  implications  of  these  data  for  both  short-pulse  engine  testing  and 
contamination  tests  to  analyze  further  the  data  in-hand.  Figure  8  shows  the  radial  profile  of 
the  nu  product  at  42.5  msec  of  the  engine  bum,  and  it  is  seen  that  the  particulate  effluent  is 
slightly  asymmetric. 
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Figure  2.  Off-axis  temporal  variation  of  particulate  na  product 
(0.1225  in.  below  centerline). 
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Figure  3.  Off-axis  temporal  variation  of  particulate  no  product 
(0.245  in.  below  centerline). 


Figure  4.  Off-axis  temporal  variation  of  particulate  no  product 
(0.3675  in.  below  centerline). 
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Figure  6.  Off-axis  temporal  variation  of  particulate  na  product 
(0.625  in.  below  centerline). 
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The  average  particle  size  parameter  (x)  was  determined  using  the  degree  of  polarization 
(P)  measurements  obtained  at  the  scattering  angle  of  90  deg.  The  angle  (<£)  between  the 
scattering  plane  and  the  polarization  vector  of  the  incident  laser  was  177  deg  (rather  than  the 
design  value  of  180  deg).  As  discussed  in  Refs.  2  and  3,  the  size  parameter  results  are 
dependent  upon  the  index  of  refraction  (77)  of  the  particle.  Since  the  particulate  material  was 
assumed  to  be  unburned  liquid  fuel  droplets,  the  index  of  refraction  used  was  appropriate 
for  that  of  liquid  species.  Figure  9  shows  the  variation  of  P  with  x  for  an  index  of  refraction 
of  1.2  for  <t>  -  177  deg.  Also  shown  in  Fig.  9  is  the  experimentally  determined  value  of  P  for 
a  typical  centerline  measurement.  It  is  seen  that,  for  this  case,  x  =  2.08  ,  which 

corresponds  to  D  =  0.45  ±  0.01/tm.  Figures  10  and  11  show  similar  results  for  77  =  1.335 
and  1.5,  respectively,  and  Fig.  12  shows  the  inaccuracy  in  x  (and  D)  caused  by  the 
uncertainty  in  77.  It  is  seen  from  Fig.  12  that,  using  tj  =  1 .335,  which  is  a  value  characteristic 
for  low  viscosity  liquids,  x  ±  Ax  =  1.70  Consequently,  the  inaccuracy  in  x  due  to  the 
uncertainty  in  77  is  on  the  order  of  ±30  percent,  which  significantly  exceeds  the  uncertainty 
in  x  caused  by  experimental  imprecision. 


Figure  9.  Variation  of  degree  of  polarization  (P)  with  size  parameter 
(x)  for  particule  refractive  index  of  1.2. 
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Size  Parameter  (S) 

Figure  10.  Variation  of  degree  of  polarization  (P)  with  size  parameter 
(x}  for  particle  refractive  index  of  1.335. 
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Figure  11.  Variation  of  degree  of  polarization  (P)  with  size  parameter 
(x)  for  particle  refractive  index  of  1 .5. 
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Using  this  procedure  for  all  the  radially  dependent  data  points  acquired  for  the  17-msec 
firing  time  enabled  the  average  diameter  D  to  be  determined,  and  the  results  are  shown  in 
Fig.  13.  From  these  results  it  is  seen  that  the  measured  size  parameter  (x)  is  on  the  order  of 
1 .75  (which  corresponds  to  a  particle  diameter  of  approximately  0.4  /im)  and  that  this  result 
is  comparatively  constant  along  the  radial  direction. 

Finally,  assuming  time-independence  of  the  particle  diameter  values  shown  in  Fig.  13, 
the  number  density  (n)  of  the  particles  was  found  from  the  measured  no  products.  Figure  14 
shows  these  results  for  each  of  the  spatially  resolved  radial  positions  as  a  function  of  engine 
combustion  time. 
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Radial  Distance,  in. 

Figure  13.  Radial  variation  of  size  parameter  fx). 

Regarding  the  uncertainty  in  the  results  shown  in  Fig.  14,  one  recalls  that  the  imprecision 
values  estimated  for  the  size  parameter  and  particle  diameter  are  less  than  ±  10  percent. 
Since  o  «  D3,  the  corresponding  imprecision  values  for  the  number  density  are  on  the  order 
of  ±  30  percent,  which  characterizes  the  relative  imprecision  ofvthe  variations  of  n  with 
respect  to  pulse  time  and  vertical  distance.  Further,  an  estimate  of  ±  100  percent  as  the  order 
of  magnitude  of  the  absolute  uncertainty  in  n  is  provided  by  recollection  of  the  uncertainty 
( ±  30  percent)  in  x  resulting  from  our  ignorance  of  the  exact  value  of  the  index  of  refraction 
(v)- 

To  obtain  an  order -of-magnitude  estimate  of  the  particulate  flux  (Jp),  one  must  compute 
the  flow  speed  at  the  exit  plane  of  the  nozzle.  For  this  purpose  it  is  assumed  that  at  the  exit 
plane  of  the  100:1  area  ratio  nozzle  the  specific  heat  ratio  (y),  the  gas  species  gram  molecular 
weight  (M),  and  the  gas  temperature  (T)  are  (Ref.  2) 

y  =  1.2, 

M  =  30  gm/mole, 

and 

T  =  87 OK 

From  these  estimates  one  finds  the  free-stream  flow  speed  to  be 

»  3x  10s  cm/sec. 
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Therefore, 


n  v 


00 


3  x  105  panicles  /cm*  /sec, 


and  the  uncertainty  in  J  is  on  the  order  of  ±  100  percent.  Since  the  assumed-liquid 
particulate  bulk  density  (eg)  is  approximately  equal  to  1  gm/cc,  one  finds  the  contaminant 
mass  flux  (Jm)  to  be 


D3/e)pg  jp 

=  10"®  gm/  cm^/  sec 


3.2  COMBUSTION  INSTABILITY 

It  has  been  shown  that  the  5-lbf  liquid  bipropellant  engine  which  was  the  subject  of  the 
study  reported  in  Refs.  1  and  2  achieved  steady  state  with  respect  to  particulate  effluent  in  a 
characteristic  time  on  the  order  of  20  to  30  msec.  However,  there  remains  the  possibility  that 
the  existence  of  combustion  instability  in  the  combustion  chamber  may  not  be  adequately 
simulated  with  a  short  pulse  duration  of  the  engine  firing.  In  particular,  if  the  characteristic 
frequencies  of  the  fundamental  acoustical  modes  of  the  combustion  chamber  are  sufficiently 
low  and  the  corresponding  characteristic  periods  are  sufficiently  large  with  respect  to  the 
engine  firing  duration,  only  a  portion  of  the  combustion  instability  will  be  observed  for 
short  pulse  firings.  Figure  15  sketches  this  behavior  for  two  cases  of  combustion  instability, 
and,  from  Fig.  15,  it  is  seen  that  inaccurate  simulation  can  result  for  a  low  oscillatory 
eigenmode  frequency  in  the  combustion  chamber  if  the  engine  firing  time  (t)  is  less  than  the 
oscillation  period.  Therefore,  it  is  worthwhile  to  evaluate  the  characteristic  frequencies  and 
periods  for  the  fundamental  oscillation  modes  of  the  combustion  chamber  and  to  compare 
these  quantities  with  the  engine  pulse  durations  to  verify  the  accuracy  of  simulation  of  short 
engine  pulses. 

To  describe  the  combustion  instability  phenomenon  the  combustion  chamber  acoustical 
model  will  be  that  of  an  open  Helmholtz  resonator.  The  characteristic  oscillation 
frequencies,  f^,  are  given  by  (Ref.  4) 


f..k  =  (ao/2)  [(a^./R)2  +  (k/L)2]'4  ,  t,  j,  t  >  0 


(1) 
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where 

ao  =  equilibrium  (isentropic)  speed  of  sound  in  the  combustion  chamber, 

R  =  combust  ion  chamber  radius, 

and 

L  =  combustion  chamber  length. 

Also,  ajj  is  defined  by  the  zeros  of  the  Bessel  function  (J„)  derivative  as  follows: 

0Jn(orR)/*l  =0, 

L  "  J  r=  R 

or 

■fn  («ij)  =  0  s  J/,  aij) 


Figure  15.  Combustion  instability  versus  engine  firing  duration. 


Equation  (1)  shows  clearly  that  the  lowest  characteristic  frequencies  are  the  fundamental 
mode  values  f  ioo,  foio.  and  foot-  If  these  three  frequencies  are  sufficiently  large  for  a  specified 
engine  pulse  duration  to  preclude  simulation  problems  caused  by  combustion  instability, 
then  the  even  larger  harmonic  frequencies  are  definitely  inconsequential.  Quantitatively,  it  is 
demanded  that 
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r  >  rijk  s  W  5  >•  j.  k  >  0 

Rewriting  Eq.  (1)  in  dimensionless  form,  one  obtains 


?ijk  -  rijk  *  <V/D>  “  [»fj  +  k2  (D/2L)2  ]"** 


(2) 


and  Fig.  16  shows  the  variation  of  for  the  fundamental  modes  as  a  function  of  the  ratio 
L/D.  Both  the  tangential  and  radial  mode  quantities  tioo  and  foio,  respectively,  are 
independent  of  L/D,  and 

roo  i  =  ^L/D 


Figure  16.  Combustion  instability  period  versus  chamber  shape  factor. 
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For  the  engine  employed  in  the  studies  of  Refs.  1  and  2,  the  combustion  chamber  diameter 
and  length  were  0.50  in.  (1.27  cm)  and  1.5  in.  (3.81  cm),  respectively.  Estimating  the 
combustion  chamber  temperature,  molecular  weight,  and  specific  heat  ratio  to  be  3000° K, 
30  gm/mole,  and  1.20,  respectively,  one  finds  that  ao  *  10s  cm/sec. 


Therefore,  one  finds  that 

r10Q  =  K  10"*  sec, 

TqIq  =  2.17  x  1CT5  sec, 

and 


r001  =  ^.6  x  10  5  sec  =  0  <10-+  sec). 

Since  the  minimum  engine  pulse  time,  Tmj„,  is  0(10"3  sec), 
then 

rmin  >:>  r001  >  r010  >  rI00> 


Thus,  no  simulation  problems  associated  with  combustion  instability  are  anticipated. 

4.0  CONCLUSIONS 

The  results  of  this  work  have  shown  that:  a)  the  mass  density  and  number  density  of  the 
particulate  effluent  of  a  5-lbf  liquid  bipropellant  engine  attain  steady-state  characteristics  in 
a  time  of  20  to  30  msec  after  initiation  of  the  opening  of  the  propellant  valve;  b)  The  particle 
average  size  was  approximately  0.4  fim  and  exhibited  radial  variation  at  the  axial  location  of 
1 . 1  in.  (2.8  cm)  relative  to  the  nozzle  exit  plane;  c)  Both  the  radial  and  temporal  variations  of 
the  particulate  number  density  (n)  were  determined,  n  was  found  to  be  on  the  order  of  1 
particle/cc,  and  a  slight  cylindrical  asymmetry  was  observed;  d)  The  particle  mass  density 
was  estimated  to  be  on  the  order  of  3  x  10“ 14  gm/cc,  and  the  order -of-magnitude  estimates 
for  mass  flux  density  (Jm)  and  number  flux  density  (Jp)  were  found  to  be  10“ 8  gm/cm2/sec 
and  105  particles/cm2/sec,  respectively;  e)  The  possibility  that  combustion  instability  might 
affect  the  accuracy  of  simulation  was  considered,  and  the  results  indicated  that  tests  using 
10-msec  engine  burn  durations  could  be  employed  with  confidence  if  the  combustion 
chamber  ratio  L/D  was  not  more  than  30. 

In  conclusion,  it  is  suggested  that  the  accuracy  of  simulation  of  pulsed  engine  testing  with 
regard  to  particle  kinetics  and  engine  performance  has  been  demonstrated. 
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NOMENCLATURE 

ao  Isentropic  sound  speed  in  combustion  chamber 

D  Combustion  chamber  diameter 

D  Mean  diameter  of  particle 

%  Helmholtz  resonator  characteristic  frequencies 

Jm  Particulate  mass  flux 

Jn  Bessel  function  of  first  kind  and  order  n 

Jp  Particulate  number  flux 

k  Acoustic  mode  wavenumber 

L  Combustion  chamber  length 

M  Gram  molecular  weight  of  gas  species 

n  Particulate  number  density 

P  Polarization  parameter 

R  Radius  of  combustion  chamber 
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r  Radial  coordinate  (cylindrical  representation) 

T  Gas  temperature 

t  Time 

v  Free-stream  flow  speed 

x  Axial  position  of  flow  field 

x  Mean  diameter  of  particle 

y  Plume  radial  coordinate  perpendicular  to  axial  centerline  of  flow 

o<ij  Zeros  of  Bessel  function  derivatives 

7  Isentropic  ratio  of  specific  heat  capacities 

ij  Refractive  index  of  particle 

X  Wavelength  of  laser 

Qt  Bulk  density  of  liquid 

Op  Particulate  mass  density 

o(x)  Mie  scattering  cross  section 

r  Engine  pulse  duration 

Tijt  Characteristic  period  =  1/fjjk 

Tjjk  Dimensionless  characteristic  period  =  (a0/D)rijj, 
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